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Abstract: In this paper, an adaptive law with an integral action is designed and 
implemented on a DC motor by employing a rotary encoder and tachometer sensors. The 
stability is proved by using the Lyapunov function. The tracking errors asymptotically 
converge to zero according to the Barbalat lemma. The tracking performance is specified 
by a reference model, the convergence rate of Lyapunov function is specified by the matrix 
Q and the control action and the state weighting are restricted by the matrix T . The 
experimental results demonstrate the effectiveness of the proposed control. The maximum 
errors of the position and velocity with the integral action are reduced from 0.4 V and 1.5 V 
to 0.2 V and 0.4 V, respectively. The adaptive control with the integral action gives 
satisfactory performance, even when it suffers from input disturbance. 

Keywords: MRAC; integral action; Lyapunov stability; asymptotical convergence 
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1. Introduction 

Various electromechanical motors have been used for industrial applications, e.g., electric motors, 
piezomotors and hydraulic actuators. Generally, it is necessary to enhance the motor performance with 
the feedback control based on state measurements. Khorrami et al, used linear motors to establish 
ultra- accurate high-speed six degree-of-freedom manipulation [1]. To measure the angle and speed of 
motor shafts, rotary encoders and tachometer sensors have been widely used because of their high 
performance and reliability [2-4]. 

As digital angle-measuring sensors, rotary encoders consist of optics, mechanics and electronics. 
Compared to analog angle-measuring sensors, digital rotary encoders have simple structures while 
preserving high accuracy. In order to achieve accurate motion control, the velocity can be measured by 
using tachometer sensors which usually consist of tachogenerators and circuits. The tachogenerator can 
also give the voltage output which is proportional to the speed of the rotational motor. 

To achieve precise motion, it is necessary to measure and input both the angle and the speed of 
motor shafts by employing a rotary encoder and tachometer sensors. Based on the angle and speed 
sensors, various controllers have been designed for motors [5,6]. Generally, motors can be controlled 
by conventional PID controllers. Nikulin and Frantsuzova presented a modified PD controller to 
ensure the desired system speed and damping vibrations [7]. Chaiya and Kaitwanidvilai provided a 
robust PID controller which could control the motor speed [8], but PID controllers have limited 
performance in the presence of disturbances and uncertainties. Xu and Yang presented a simple and 
robust speed control scheme for a permanent magnet synchronous motor to enhance the performance 
robustness [9]. Nouri et al, proposed a model-following adaptive controller for the speed control of a 
motor drive system [10]. Melkote and Khorrami proposed adaptive control for direct drive brushless 
DC motors [11]. Moreover, intelligent algorithms have also been investigated [12-15]. Fallahi and 
Azadi added neural network sliding mode control to enhance the adaptive control of motors [16]. 
However, intelligent controllers are usually very complex. It is difficult for engineers to design and 
optimize intelligent controllers, e.g., Xu and Huang [17] designed an iterative learning controller 
(ILC), but found that the reference signal had to be pre-filtered in order to satisfy the complex initial 
conditions of ILC. Furthermore, several accurate models were identified at different voltage ranges, 
and the iterations were implemented offline. 

This paper uses Model Reference Adaptive Control (MRAC), which is a convenient approach to 
satisfy the requirements of designers [18]. The general idea of MRAC is to create a closed loop 
regulator with parameters that could be updated to match a desired response. The desired performance 
is specified by a stable reference model, and the parameters of the adaptive law are adjusted based on 
the errors between the reference model and the plant, as shown in Figure 1 . 

Though the adaptive control has shown its effectiveness in achieving robust performance without 
the knowledge of parameter values, a comprehensive design approach is still necessary for engineering 
applications. In order to obtain high-performance adaptive control in the presence of disturbances, this 
paper presents a comprehensive design approach in which both the bandwidth and damping ratio can 
be included in the proposed controller. 
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Figure 1. MRAC control sketch of the motor. 
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In this paper, two kinds of MRAC are designed and then used on a motor by employing a rotary 
encoder and tachometer sensors. The tracking error can converge to zero with the integral action in the 
presence of input disturbances. The experimental results are presented to investigate the effectiveness 
of the proposed control approach. 

2. Adaptive Control Design without the Integral Action 

2.1. Angle and Speed Sensors 



Figure 2 shows the experimental setup. It consists of a DC motor, a computer with LabView 
software, a drive interface, an amplifier, a tachometer sensor and a rotary encoder. The tachometer 
sensor is used to measure the rotation speed of the motor shaft. The output voltage of the tachogenertor 
{i.e., tachometer generator — a device to measure the rotational rate of the motor shaft with the 
internally generated electrical signal generated by the motor shaft) is proportional to the motor speed. 
Then, the voltage is applied to the voltmeter in which the dial can be calibrated in speed units 
{i.e., usually in revolutions per minute (rpm)). In addition, a rotary encoder is used to measure the 
linear rotation angle of the motor draft. The rotary encoder can convert the ration angle to digital 
voltage. Thus, both the angle and speed can be measured and controlled. 



Figure 2. Sketch of the motor experiment. 



Sensors 2013, 13 



4745 



After calibration, the generated voltage is a factor of the motor speed. Then, the motor speed can be 
fed to the adaptive controller. Next, Figure 3 is given to present the working principle of the rotary 
encoder in which phase A, phase B, and phase Z are the three-phase models to represent the output 
signal counts for an incremental encoder [19]. The rotary encoder consists of the main grating, index 
grating, illuminant and photosensitive device. In this paper, an incremental encoder is used because of 
its simplification, small size and high-speed response. There is one impulse corresponding to every 
grating. The summation of impulses represents the angle position of the motor. 

Figure 3. Working principle of the rotary encoder. 
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2.2. Adaptive Controller Design and Stability Proof 

This section presents a design method of the adaptive control without the integral action. The 
stability proof is also presented. To achieve good tracking performance, a MRAC is designed to drive 
the tracking error to zero. Considering simplification, the transfer function in Equation (1) is used to 
give the motor dynamics: 

K 

(1) 



G = 



where K denotes the DC gain of the velocity transfer function, and r denotes the time constant, r are 
unknown parameters, but their signs can be tested in the experiment. The state space model of 
Equation (1) can be rewritten in Equation (2): 

x = A p x + gbu 

y -Cx 



(2) 



where: 



0 1 



0 
1 



,g = K, C 



1 0 
0 1 



To specify the desired performance, this paper employs a stable reference model as shown in 
Equation (3), from which the domain and frequency index can be specified: 
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GL =■ 



s 2 + 2^<y n 5 + col 



(3) 



where £ is the damping ratio and co n is the natural frequency of the reference model. The corresponding 
state space description can be given as: 



(4) 
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Equation (4) can be further written as: 



K = K x m + H m br 



The following non-adaptive control law is used: 

u = 07x p +6*r 

where 0* x and 0* are the exact gains of the controller Equation (5). 
Then, the closed loop system can be given by: 

x p = [A p + gbe; T }c p +(g0:)br 



(5) 



(6) 



(7) 



The feedback control system achieves the performance as the matching condition below: 

A P + gb&7 = A„ 

The exact control gains 6*J and 9* r guarantee that the closed loop system matches the reference 
model. Actually, the exact values of 6*J and 0* r are unknown. The controller in Equation (5) can be 



rewritten to: 



u = e T (x)x+0r{t) 



where 0 T X and 0 r can be determined by the adaptation law in Equation (9): 

0 = 0* + 0 = 20 = -2sgi( g)rxe T Pb 

where 0* is a constant or the slow time-variant parameter (i.e., 0* =0). 
The adaptive control law is nonlinear, as shown in Figure 4. 



(8) 



(9) 



Figure 4. Adaptive control block. 
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To give the tracking performance, the parameter errors <j) x and <f) r are defined by: 
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<l> r =e r -e; 

Then, the closed loop system can be rewritten to: 

*p = A p x p + gb(0 T x (x)x p + 0 r r{t)) 
= A m x p + gbf x x p +gb0 r r(t) 

Next, let e be the state error and e-x p -x m . 



(10) 



(11) 



By comparing the closed loop system in Equation (1 1) with the reference model in Equation (5), the 
dynamics of the tracking error e can be given: 



e = A m e + gbf x x p + gb9 r r(t) 

= K e + sH T x 



(12) 



where: 







~ X P~ 




,x = 




A 




r 



0= 



In the reference model in Equation (4), a stable matrix is used for A m that can satisfy the algebraic 
Riccati Equation (13) (i.e., for any symmetric positive definite matrix Q, there is a symmetric positive 
definite matrix P satisfying Equation (13)): 



To prove the stability, a Lyapunov function candidate is used as follows: 

V(e,0) = e T Pe + \g\fr- l <f> 



(13) 



(14) 



where r is a positive definite matrix. 

The time derivative of the Lyapunov function candidate can be given: 

V(e,0) = 2e T Pe + 2gfr- 1 <j> 

= e T {A T m P + PAJe + 2ge T Pbf + IgfT' 1 ^ 
= -e T Qe + 2ge T Pbf + 2gf r~ l ^> 

Then, Equation (15) can be rewritten to: 

V(e,<p) = -e T Qe + 2ge T Pb<fi T x-2sgi(g)\g\fxe T Pb 
- -e T Qe < 0 



(15) 



(16) 



where g =sgn(|g|)|g|, <j) T x is a scalar. 

Thus, the system is stable according to the Lyapunov theorem [20]. Moreover, ||e||, e,||^||, x n , 0 X 
and 0, are bounded, and: 



\e T Qedt = V 0 -V x <V 0 <c x 
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According to the Barbalat lemma [20] here V : [0,+qo] R is a uniformly continuous function on 
[0,+oo] . Supposing limJV \z)dr exist and be finite. Then, there is V — > 0. 

f-»oo() t— >KJ 



Here e is bounded. Furthermore: 



According to the Barbalat lemma, the tracking error e is asymptotically stable and lime = 0 
(i.e., lim(||x m -x„||) = 0). 

To improve the adaptation rate, the proportional action can be added to the adaptive law [14]. 
Finally, the adaptive law can be given in Equation (17): 

6 = -(ft | sgn(£ p )Tae x dt + y 2 sgn(£ p )Tae x (17) 
where y x = 1 and y 2 = 0.5 , and they are used in experiments. 

2.3. Experimental Studies of Adaptive Control without the Integral Action 

Experimental studies of the adaptive control are now presented. Firstly, p can be used to give the 
convergence rate of the adaptive law, as shown in Equation (18): 



V 



e T Qe 



_^e T Qe ^_{Q) 
e T Pe + \g\<?> T r- l <f> e T Pe A mn (P) 



(18) 



After several trials, the reference model is given: 

0 f 



K = 



■25 



, o) = 5, £ = 0.9 



and the matrix Q is selected as: 



Q = 



10 0 
0 1 



The selection of Q means that the position tracking is more important than the velocity tracking. 
Moreover, the sampling interval h should be less than the max interval in order to handle the fastest 
adaptation rate [18]: 

1 1 _ 1 ^ n (P) 

20 p m 20 A^Cfi) ( 19 ) 



h< 



The matrix P is solved by the ARE in Equation (13): 

"3.74 0.2 



P = 



and denoting r 



r = 



0.2 0.078 

1 0 0^ 
0 0.5 0 
0 0 1 
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The tracking performance of the adaptive control is shown in Figures 5 and 6. It can be found that 
the tracking errors of angular position and velocity asymptotically converge to zero in one period, and 
the maximum tracking errors of position and velocity are 0.4 and 1.5, respectively. 

Figure 5. The reference position, the motor position and the tracking error (i.e., (top), 
(middle) and (bottom), respectively). 




Figure 6. The reference velocity, the motor velocity and the tracking error (i.e., (top), 
(middle) and (bottom), respectively). 




To guarantee the stability of the adaptive control, it is necessary to demonstrate the boundness of 9 
Figure 7 shows the responses of 6 . It can be seen that the estimation of 9 = [ 6^ 6 Xi 6, ] is bounded. 
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Figure 7. 9 responses under square wave tracking (top) and sinusoidal tracking 
(bottom), respectively. 




2.4. Investigations of Q andT 

Q and r are two important parameters of the adaptive controller. Different values of Q and r are 
adopted to investigate their influences, here two group values of them are: 



Ql = 

a = 



10 0 
0 1 
30 0 

0 30 

1 0 0 
0 0.5 0 
0 0 1 



3.74 0.2 
0.2 0.078 
48.7 0.6 
0.7 1.778 
10 0 0 



r 2 = 



5 0 
0 10 



where Q 2 is larger than Q l . 

The adaptation rates can be respectively estimated: 



^max 1 



P 2 = 



l mn (P)_ 0.0669 
^«(fi)_ 30 



= 149.5 



^ a (P) 1-77 



= 16.9 



where p is the adaptive rate, and p l is faster than p 2 . 

During experiments, the adaptive control using Q ] and Q 2 have similar tracking errors, but the 
control signals are different. Figure 8 shows the voltages of the adaptive control. The control voltage 
of the adaptive control with (Q 2 , r 2 ) is four time larger than that of the adaptive control with (Qj, /";), 
though there is no significant differences in the tracking errors. 
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3. Adaptive Control with the Integral Action 

3.1. Adaptive Controller Design with the Integral Action and Stability Validation 

To further investigate the adaptive control in the presence of disturbances, this section presents the 
adaptive control with integral action. The integral control action can be added to the adaptive law 
through Equation (20): 

Xj=x x -r (20) 
where r is the reference signal and x, is a new state: 



Xj = - rdz 



(21) 



The augmented plant can be rewritten to: 

"0 1 0 



(x^ 



yX^j 



0-0 
r 

1 0 0 
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"0" 




"0 




+ g 
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u + 


0 






0 




-1 



(22) 



Similar to Section 2.2, the control law can be given by: 



(23) 



Then the closed loop system can be given as: 
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(24) 
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where: 
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0 




a 22 
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0 


0 



i p + gb0 x gb9 I 



0 



According to the matching condition in Equation (7), the following equation can be given: 



(25) 



A p + gb9 x gb9] 
1 0 



(26) 



The reference model can be given: 



X m ~ An X m + 



0 



1 



(27) 
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Equation (29) can be rewritten to: 

X p = Ar, X p + 

The error equation can be given 



gb(f> x gbfa 
0 0 





"0 " 


X P + 






-1 


p 





Choosing the Lyapunov function as Equation (32) 



V = e 1 Pe + \g\p 1 r-'p 
Then, the time derivative of V can be obtained: 

^ = 2e T Pe + 2\g\p T r- l p = e T (A*P + PAJe + 2ge T Pbif> T x p + 2\g\y T T^y 



(28) 



Actually, the exact value of 0 X and 0\ cannot be known, so the corresponding estimated values of 
G x and 0, are used instead. We note that (j> x and fa are the estimate errors as shown in Equation (28): 

fa =9,(1)- 9] 

and 



and e - x - x„ 



where 9 x (t) and 0j(t) are the estimated values 9* x and 0] , respectively. 

Substituting u into the plant model in Equation (22), the closed loop can be given: 

"0 



(29) 



(30) 



(31) 



(32) 



(33) 



Denote that the adaptive law is: 
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4(0 



= -sgn(g),T 



(34) 



Then, P can be gotten from the ARE Equation (35): 

^P+PAn=~Q (35) 

where Q is a positive definite symmetric matrix, so the following Equation (36) can be obtained: 

dV 



dt 



= -e'Qe<0 



(36) 



It can be concluded that the closed loop is Lyapunov stable, and e and <p are bounded. As in 
Section 2.2, the tracking error e is asymptotically stable. 
The adaptive control law is shown as follows: 



u = 0*x p + 0i*i 



6 = 



6, 



= -sgn(g)r 



x.. 



Pb 



(37) 
(38) 



where <j> « 6 , P can be computed through the ARE. 

3.2. Reference Model Design 

This section presents how to specify the requirements through the reference model G m . Firstly, the 
third order model in Equation (39) is given to satisfy the matching condition of the adaptive integral 
control. To specify performance indexes such as the bandwidth and damping ratio, the reference model 
G m can be approximated by the second order model G ml : 



f~< f f 

s +2£w n +w n s + a„ 



a m w 2 n 



(39) 



(s + P^s + P^s + aJ 

where ||a m || > SjjpJ (i.e., p l and p l are dominant poles which mainly contribute to the dynamic 

response of the reference model.) 

To satisfy the matching conditions, the reference model is: 

i- 0 
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x 2 
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K X 3J 
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0 


\ X 3J 
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y = [l 0 0] 



(40) 



Both Equations (39) and (40) describe the same reference model. The state space Equation (40) is 
transformed into the transfer function: 
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D + C(sI-A) l B = [l 0 0] 
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(41) 



23 



Compared with the transfer function in Equation (39), it is seen that: 



(42) 



a 23 = -(y„a,„ 



The dynamics of the reference model should match the dynamics (i.e., natural frequency) of the 
motor system and the sampling capability. After testing several groups of values, C, = 0.9 , w n = 3 and 
a m = 9 are chosen in this paper. Then, there are a 2l = -57.6, a 22 = -14.4 and a 23 = -81 . The reference 
model can be rewritten to: 



V 




x 2 









0 1 
-57.6 -14.4 

1 0 



0 
81 
0 



y = [l 0 0] 



(43) 



The step responses of G ml and G m are shown in Figure 9. It is seen that G ml contributes most of the 
response for the third order reference model G m . The fast dynamics of the reference model is 
influenced by the part aJis + a m ): 

Figure 9. Step responses of G m (solid line) and G ml (star line). 
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3.3. Experimental Studies on the Adaptive Control with Integral Action 



In Section 2, it has been found that a large Q results in fast convergence while needing strong 
action and fast sampling rate. In contrast, a small Q results in slow convergence and cannot 
give better performance. A suitable Q should be given to satisfy the convergence rate and 
hardware limitation. 



After trials, Q is adopted: 



Then, P can be computed: 



In addition, r is adopted: 



Q = 



15 

0 

0 



0 0 
10 0 
0 15 



P = 



32.9913 0.7111 34.4582 
0.7111 0.3966 0.0926 
33.4582 0.0926 62.9307 



r = 



o o 

2.5 0 
0 4 



The experimental results of the position and velocity tracking are shown in Figures 10 and 11, 
respectively. The maximum tracking errors of position and rate are less than 0.2 V and 0.4 V, respectively. 



Figure 10. The reference position, the motor position and the tracking error {i.e., (top), 
(middle) and (bottom), respectively). 
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Figure 11. The reference velocity, the motor velocity and the tracking error (i.e., (top), 
(middle) and (bottom), respectively). 
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As shown in Figure 12, the magnitude of the controller is less than 0.85 V, and there is no chatter 
phenomenon. Compared to the former adaptive control without the integral action, the case with the 
integral adaptive control can reduce the reference errors of position and velocity by as much as 50%. The 
experimental results demonstrate that the integral action is effective to improve the tracking performance. 



Figure 12. Control signal. 




3.4. Tracking Performance with Input Disturbance 

In order to investigate the performance of the adaptive control, the input disturbance is considered 
in this section. The square wave disturbance with magnitude 1 is adopted. The reference is also square 
wave. During the experiment, the disturbance is added by using a LabView block as an input 
disturbance. The disturbance enters the close loop systems as the reference signal enters. The tracking 
errors of position and velocity are shown in Figure 13. It is seen that the maximum position and 



Sensors 2013, 13 



4757 



velocity tracking errors are less than 0.2 and 1, respectively. Also, it can be found that there is no 
difference between the two results of the position tracking errors shown in Figures 10 and 13, 
respectively. The position error is the most important one for the motor tracking in this paper, thus it 
can be concluded that the degradation of tracking performance is not significant. The experimental 
result indicates that the adaptive control with the integral action suppresses disturbance and tracks the 
reference simultaneously. 



Figure 13. Tracking errors of position (top) and velocity (bottom). 
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3.5. Comparison between Experiment and Simulation Results 

Finally, a comparison between experiment and simulation results is presented in this paper. 
According to the motor manual, the time constant is 0.25. The DC gain is 5.2. Figure 14 shows the 
comparison between experiment and simulation results. 

Figure 14. Comparison of experiment and simulation results (case I: experimental results 
of Adaptive control; case H: experimental results of adaptive control with integral action; 
and case III: Simulation results of simulation of adaptive control with integral action). 
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The simulation and experiment have similar responses. The tracking errors of the position and 
velocity from experiment study are larger than that from simulation study. In the three cases, all the 
tracking errors asymptotically converge (i.e., the stability can hold). The simulation of the proposed 
adaptive control with the integral action gives the best performance. The tracking error from 
simulation study can converge to zero, but the tracking errors from experiment study cannot converge 
to zero due to the friction, the sensor noise and the motor dead-zone. In the experiment study, the 
performance degradation is acceptable. 

4. Conclusions 

In order to accurately control motors, this paper employs both a rotary encoder and tachometer 
sensors to measure the angle -position and speed, respectively. Based on the measurements, two 
adaptive controllers are developed for the motor system. The stability and convergence are validated 
by the Lyapnov theorem and Barbalat lemma. Then, the control system is implemented by using 
Labview. Experimental results indicate that the tracking errors of the motor position and velocity 
asymptotically converge to their error ranges. In the presence of disturbances, the adaptive controller 
with integral action presents better performance than the case without integral action. 

State estimation approaches are encouraged to investigate for possible output feedback control. In 
this paper, the proposed adaptive controllers needs full state information, but actually some states (e.g., 
motor velocity) are not provided. Thus, the state estimation must be designed in the possible output 
feedback control. 
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